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ABSTRACT 


r 

A  general  method  for  solving  the  inverse  diffraction  problem  is 

presented.  It  is  based  on  an  identity  of  Bojarski  which  states  that  y  (x) 

and  r(p)  are  a  Fourier  transform  pair.  Here  y(x)  is  the  characteristic 

function  of  the  target  (y  =  1  inside  the  target,  y  =  0  outside), 

2w 

p  =  —  ^  1  is  the  frequency,  ^  is  a  unit  vector  specifying  the  aspect, 
”  c 

and  r(p)  can  be  obtained  by  measurement  of  the  backscattered  electromag¬ 
netic  far-field  at  frequency  =  X  1p|  and  aspect  ^  =  jpj  p  .  If  data 
is  obtained  in  any  subset  D  of  p-space,  the  method  yields  partial  or  com¬ 
plete  information  about  the  target  geometry.  It  is  used  to  rederive  earlier 
results  very  simply,  and  to  obtain  a  significant  new  solution,  in  which  the 
target  geometry  is  completely  determined  using  frequencies  only  in  a 
practical  frequency  band  and  aspects  in  a  narrow  cone.  , 
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SECTION  I 


INTRODUCTION 


The  inverse  diffraction  problem  requires  the  determination  of  the  size 
and  shape  of  a  scattering  target  from  given  data.  For  example,  the  far- 
field  backscattered  at  certain  aspects  and  frequencies  may  be  given.  Any 
approach  to  the  problem  must  be  based  on  either  the  exact  or  some  approxi¬ 
mate  theory  of  direct  scattering. 


As  in  other  approaches  to  the  inverse  problem,  the  theory  developed 


here  will  be  based  on  the  physical  optics  or  Kirchhoff  approximation 


[1] 


for 


direct  scattering.  Our  starting  point  will  be  a  remarkable  identity  recently 

[2I 

obtained  by  N.  Bojarski.  A  derivation  of  the  identity  is  given  in  Section 


II.  There  we  introduce  a  three-dimensional  space  of  vectors  p  whose 
direction  coincides  with  the  aspect  direction  and  whose  magnitude  is 
p  =  |p|  =  _  where  co  is  the  frequency.  We  also  introduce  the  charac¬ 
teristic  function  y  (x)  of  the  target  and  a  function  p(p)  which  can  be  obtained 
by  measuring  the  far-field  backscattered  in  the  direction  of  p  at  the  fre- 
quency  co  =  -  Ipj  .  The  function  •y(x)  is  1  inside  the  target  and  0  outside 
of  it.  Then  Bojarski 's  identity  states  that  the  functions  y(x)  and 
r(p)  =  (p^)  ^  ^p(p)  +  p  (-p)J  are  a  Fourier  transform  pair.  Thus 

if  p(p)  could  be  measured  in  all  of  p- space  we  could  immediately  obtain 


y  (x)  and  hence  the  target. 


In  practice  p(p)  and  hence  r(p)  can  be  measured  only  in  a  restricted 
range  of  aspects  and  frequencies,  i.e. ,  in  a  subset  D  of  p-space.  In  the 
general  theory  presented  in  Section  in  we  show  how  these  restricted  scat¬ 
tering  data  can  be  used  to  obtain  information  about  the  geometry  of  the 
target.  In  some  cases,  even  when  D  is  a  relatively  small  subset  of  p-space. 


the  target  can  be  completely  determined.  In  other  cases,  partial  but  often 
useful  information,  such  as  the  width  function,  cross-sectional  area  function, 
volume,  etc. ,  can  be  obtained. 

This  suggests  a  possible  reformulation  of  the  inverse  diffraction  prob¬ 
lem  itself.  In  the  usual  formulations  some  data  domain  D  is  prescribed 
and  one  asks  for  proofs  of  the  existence  and  uniqueness  of  a  target  that 
produces  the  specified  data,  as  well  as  for  constructive  algorithms  which 
yield  the  target  geometry.  In  the  reformulated  version,  for  each  prescribed 
data  domain  D  ,  one  asks  for  the  maximum  determination  of  geometric 
parameters  that  can  be  obtained  and  for  algorithms  for  constructing  these 
parameters.  For  relatively  "small"  domains  D  these  parameters  may  not 
determine  the  target  uniquely.  In  such  cases  we  obtain  only  the  class  of 
targets  (e.  g. ,  all  targets  with  a  given  volume)  with  the  specified  parameters. 
For  unique  determination  of  the  target,  or  for  a  given  type  of  partial  speci¬ 
fication,  one  may  also  ask  for  the  "minimum"  or  most  "accessible"  data 
domain  D  needed. 

In  Section  IV  we  apply  the  general  theory  to  four  elementary  examples. 

In  these  examples  D  is  first  all  of  p-space,  then  a  plane,  a  line,  and 

finally  a  point  in  p-space.  To  readers  familiar  with  physical  optics  methods 

for  the  inverse  diffraction  problem  the  nature  of  the  results  will  be  familiar, 

but  the  ease  and  simplicity  with  which  they  are  obtained  demonstrate  the 

power  of  the  present  method.  In  Section  V  we  obtain  a  complete  target 

specification  using  data  in  and  infinitesimally  near  an  infinite  plane.  Since 

this  plane  passes  through  the  origin,  it  contains  all  frequencies,  and  all 

aspects  in  some  "equatorial  plane"  of  the  target.  Here,  too,  the  results  have 

r  3l 

been  obtained  earlier  by  a  different  method  but  are  now  derived  with  great 
simplicity  and  elegance. 
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Existing  radar  techniques  are  restricted  to  a  frequency  band  of 
moderate  width.  This  band  corresponds  in  p-space  to  the  annular  region 
between  two  concentric  spheres  of  minimum  and  maximum  frequency. 

Even  if  this  annular  region  could  be  greatly  enlarged  by  future  technologi¬ 
cal  developments,  there  would  remain  a  "forbidden"  low-frequency  sphere 
about  the  origin  in  which  the  physical  optics  approximation,  which  forms 
the  basis  of  our  theory,  fails.  It  is  clear  that  to  obtain  a  practical  method, 
the  domain  D  must  lie  within  the  annular  region.  Those  examples  in 
Sections  IV  and  V  that  yield  useful  target  information  all  violate  this  require¬ 
ment.  Therefore  in  Section  VI  we  present  an  application  of  the  general 
theory  which  yields  a  complete  determination  of  the  target,  in  which  the 
domain  D  lies  entirely  within  the  annular  region.  It  consists  of  a  disc¬ 
shaped  portion  of  a  plane  which  is  tangent  to  the  inner  sphere  and  truncated 
by  the  outer  sphere.  Thus  the  frequencies  lie  within  the  usable  band  and 
the  aspects  lie  within  a  cone  of  semi-angle  6  .  For  a  sufficiently  large 
upper  frequency  limit,  6  can  be  made  arbitrarily  small.  Thus  the  results 
of  Section  VI  might  form  the  theoretical  basis  for  a  workable  radar 
recognition  system. 

We  note  that  since  r(p)  =  2\/ir{p^)  ^  [p(p)  +  p  (-p)J  ,  measure¬ 
ment  of  r(p)  in  a  domain  D  requires  measurement  of  p  in  D  and  in  D'  , 
which  is  the  reflection  of  D  through  the  origin.  Thus  if  D  corresponds  to 
aspects  near  the  "front"  of  the  target,  p  must  unfortimately  be  measured 
also  near  the  "back.  "  Since  this  would,  in  many  applications,  be  a  severe 
limitation,  we  have  attempted  to  eliminate  the  domain  D'  .  In  Section  VII 
we  present  an  extension  of  our  general  theory  which,  under  suitable  con¬ 
ditions,  yields  the  "front  half"  of  the  target  using  measurements  of  p  in  D. 
This  extension  can  be  utilized  in  connection  with  the  method  of  Section  VI. 

In  Section  VIE  we  derive  a  simplified  version  of  our  general  theory  suitable 


3 


for  targets  which  are  bodies  of  revolution.  By  introducing  cylindrical  coor¬ 
dinates  we  can,  in  effect,  reduce  the  dimensions  of  both^-space  and  x~space 
from  three  to  two.  Some  elementary  examples  of  this  simplified  theory  are 
presented. 

Although  we  have  used  the  general  theory  to  obtain  old  results  easily, 
and  have  obtained  one  significantly  new  result  in  Section  VI,  we  have  not 
begun  to  exhaust  the  potential  usefulness  of  the  method.  There  is  good  reason 
to  hope  that  other  useful  applications  of  the  theory  and/or  generalizations  of 
the  method,  either  by  the  present  author  or  by  others,  may  be  given  in  sequels 
to  this  paper. 
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SECTION  II 


BOJARSKI'S  IDENTITY 

The  identity  which  forms  the  basis  of  our  general  theory  is  derived” 
from  a  representation  for  the  electromagnetic  far-field  backscattered  from 
a  perfectly  conducting  target.  We  consider  an  incident-plane,  time- 
harmonic,  electromagnetic  wave  with  electric  vector. 


E .  (t  ,  x)  =  E 
— 1  -  — o 


^ik(I  •  X  +  R)  -ioJt 


Here  co  is  the  frequency  and  k  =  co/c  is  the  wave  number.  The  constant 
R  is  the  range,  I  is  a  unit  vector  in  the  direction  of  incidence,  and 
is  a  constant  vector  orthogonal  to  I  .  According  to  the  Kirchhoff,  or 
physical  optics  approximation  [  see  (3. 41)  of  Reference  l]  ,  the  backscattered 
far-field  at  the  point  x  =  -RI  is  given  by 


where 


E 

—  s 


2ikR-iaJt 


(2) 


-ik  C  ^  2ikl  •  X  , 

p  =  —  \  I  •  n  e  -  -  dS  (x)  . 

n/t 

J_i 


The  surface  integral  in  Equation  (3)  is  over  the  illuminated  portion  L  of 
the  target,  and  n  is  the  outward  unit  vector  normal  to  the  surface.  The 
power  cross  section  a  is  related  to  p  by  the  equation. 


a  =  477  R^  1e  J^/Ie  1^  =  1p1^  .  (4) 

S  1 
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We  introduce  the  vector 


p  =  -21d  (k  =  ^  ^  oj  (5) 

which  is  defined  for  non-negative  frequencies  a; .  The  magnitude  of  p  is 
p  =  2k  =  2  (u/c,  and  its  direction  is  that  of  the  unit  vector  J  - 1 
directed  from  an  origin  in  the  target  to  the  transmitter-receiver.  From  (3) 
we  see  that 


P  =  P(P)  = 


-ip  •  X  ^ 

p  •  n  e  -  -  dS(x)  , 


p  •  n 


Here  we  have  replaced  the  restriction  to  the  illuminated  region  L  by  the 
condition  p  •  n  >  0  .  The  replacement  is  strictly  valid  for  all  p  only 
if  the  target  is  convex.  In  general  we  introduce  an  additional  approximation 
in  replacing  (3)  by  (6). 

We  now  note  that 


P(P)  +  P*  (-P)  =  \  P  •  n  e”^?  -dS  .  (7) 

2-Jn 

Here  the  star  denotes  the  complex  conjugate,  and  the  integral  is  over  the 
entire  surface  S  of  the  target  B.  It  follows  from  the  divergence  theorem 
that 
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y(x)  =  (27r)  ^  j  r(p)e'?‘  -dx  .  (11) 

If  the  back-scattered  field  could  be  measured  at  all  frequencies  and  all 
aspects  _I  ,  then  p(p)  would  be  known  for  all  p  ,  and  (11)  would  yield  the 
full  solution  of  the  inverse  diffraction  problem,  i.  e. ,  the  size  and  shape 
of  B.  In  general  p(p)  is  measurable  only  for  restricted  frequencies  and 
aspects,  i.  e, ,  in  a  restricted  portion  of  p-space.  The  general  theory 
presented  in  the  next  section,  and  illustrated  in  succeeding  sections,  shows 
how  the  restricted  information  about  p(p)  can  be  used  to  yield  complete  or 
partial  information  about  the  target. 


SECTION  III 

THE  GENERAL  THEORY 

Suppose  r(p)  can  be  obtained  [  by  measurement  of  p(p)]  in  some 
portion  D  of  p-space  .  A  convenient  way  of  stating  this  is  that  we  can 
measure  the  function  K(p)  r(p)  where  K(p)  is  the  characteristic  function 
of  D  (i.e. ,  K  1  in  D  ,  K  =  0  outside  D).  Alternatively,  K(p)  could 
be  any  function  which  is  non-zero  in  D  and  zero  outside  D  .  We  introduce 
the  transform  pair 

K(x)  =  (27r)'^  Jk(p)  e^?  ‘  -dp  ,  K(p)  =  J  ‘  -dx  . 

(12) 

Since  K(p)  r(p)  can  be  measured,  we  can  construct  the  function 
f(x)  =  (27r)'^|K(p)  r  (p)  e'B  ‘  ?dp 


-2  -5/2  f  ^(P)  r  *  T  ip  .  X 

=  2  V  J  J^p(p)  +  p  (-p)jet:  ?dp  .  (13) 

P 

On  the  other  hand,  it  follows  from  the  convolution  theorem  for  Fourier 
transforms  that 


=  ^'y(x')K(x  -  x')  dx'  =  ^K(x')y(x  -  x')  dx' 


Since  f  can  be  constructed  from  measurement  and  3C  is  known  from  (12), 
(14)  is  an  equation  for  y  (x)  .  We  shall  see  that,  in  some  cases  it  can  be 
solved  to  determine  y  (x)  ,  i.e.,  the  target  B.  In  other  cases  it  yields 
partial  information  about  B.  It  is  not  surprising  that  (14)  can  sometimes 
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be  fully  solved  even  when  K  vanishes  in  much  of  p-space,  because  we  know 
in  advance  that  y  is  a  characteristic  function.  In  fact,  since  y  is  com¬ 
pletely  determined  by  a  two-dimensional  manifold,  i.  e.  ,  the  surface  of  B  , 
we  might  expect  that  (14)  can  be  solved  even  when  r(p)  can  be  measured 
only  on  a  surface,  i.  e. ,  a  two-dimensional  manifold  in  p-space.  This 
expectation  will  be  confirmed  in  Sections  V  and  VI. 

It  is  interesting  to  note  that  if 

K(-p)  =  K*  (p)  (15) 


then  K  (x)  is  real  and 


+ 


*  -ip  • 
P  (P)e  - 


K(p)p(p)  ip  .  ^ 

— ■ 

P 

Then  the  determination  of  K  and  f  is  somewhat  simplified. 


(16) 


In  general  we  proceed  by  choosing  a  domain  D  in  p-space  and  a 
convenient  function  K(p)  which  is  non-zero  only  in  D.  Then  we  determine 
K  from  (12)  and  f  (in  terms  of  p)  from  (13).  Finally,  we  use  (14)  and  the 
fact  that  y  is  a  characteristic  fimction  to  obtain  a  partial  or  complete 
determination  of  B.  The  procedure  will  be  illustrated  in  Sections  IV,  V, 
and  VI.  It  should  be  emphasized  that  those  cases  that  yield  only  a  partial 
determination  of  B  should  not  be  considered  a  failure  to  solve  the  inverse 
diffraction  problem.  Rather,  it  might  be  considered  remarkable  that  partial, 
and  often  useful,  information  can  be  obtained  from  inadequate  data. 
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SECTION  IV 


ELEMENTARY  EXAMPLES 

In  this  section  we  will  apply  the  method  of  Section  in  to  several  simple 
cases.  The  results  are  not  of  practical  importance,  primarily  because  the 
domain  D  includes  low-frequency  regions  where  physical  optics  fails  and 
both  high  and  low-frequency  regions  where  p  cannot  be  measured  by  exist¬ 
ing  techniques.  Nevertheless,  the  examples  are  useful  in  clarifying  the 
application  of  the  theory.  In  each  case  we  give  first  the  domain  D  and 
K(p).  Then  we  determine  X(x)  from  (12)  and  f(x)  in  terms  of  p  from 
(13)  or  (16)  when  (15)  is  satisfied.  Finally  we  evaluate  the  right  side  of  (14) 
and  discuss  the  resulting  information  about  the  target  B. 

Example  1: 


D: 

all  of  p- space;  K(p)  =  1;  K(x)  -  6(x)  , 

(17) 

f(x) 

-1  -5/2  r  ip  '  X  , 

=  2  TT  (te  \  e  -  -  dp  , 

(18) 

p 

f(x) 

=  ^y(x')  6(x  -  x')  dx'  =  y(x)  . 

(19) 

Here  y  (x)  is  completely  determined  by  eliminating  f  in  (19)  and  (18). 
Example  2: 

D:  the  plane  Pj^  =  0 ;  K(p)  =  <5(p^)  (20) 

K(x)  =  (27r)'M  e  ^  ^  dp^dpg  =  (27r)  "  6(X2)  a(x3)  .  (21) 
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'(o,  P2  .  P3)  ^3^3^ 


2  2 

P2  P3 


dp2dp3  , 


(27r )  ^  j'y  (x')  6  (x^  -  x^')  6(x^  -  x^')  dx' 

(27:)'^  jy(x^'  ,  x^  ,  Xg)dx^'  =  (277)"^  w  (x^  ,  x^). 


Hgts  w  ^^2  ’  ^3^  denotes  the  "width”  of  the  target  at  the  point  ^^2  ’  ^3^  ’ 
i.  e. ,  the  length  of  that  part  of  the  line  parallel  to  the  x^-axis  and  passing 
through  the  point  (o,  x  ,  x  )  which  is  contained  in  the  target.  Thus  a 
partial  description  of  B  is  given  by  (23)  and  (22).  The  description  is  incom¬ 
plete  because  many  different  targets  can  have  the  same  width-function  w. 


Example  3: 


D;  the  line  P2  ~  P3  ~  ~  ^(^2)  ^(^3) 

-3  r  -2  /  \ 

K(x)  =  (277)  J  e  dp^  =  (277)  6(x^)  , 


-1  -5/2  ^ 
f(x)  =  2  77  tfle 


p{p^>  0  ,  0)  ip^x^ 


■2 -  ^  ’ 


=  (277  )"^  j'r(x')  6(x^  - 


=  (277  )'^  Jy  (x^  ,  Xg'  ,  Xg’)  dX2’  dx. 
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Here  denotes  the  cross-sectional  area  of  the  target,  i.  e. ,  the  area 

of  that  part  of  the  plane  =  constant  which  is  contained  in  the  target. 
Thus  (26)  and  (27)  yield  a  partial  description  of  B. 


Example  4; 


D;  the  point  p  =  (Pi  ’  P2  ’  P3)  "  (^0  ’  ^  ’ 

K(P)  =  «(P2)H‘'3)^(Pi  -  Po)’ 


K  (x)  =  (27r)  e 


-3 


-2  -5/2  P(P0  ^  P  (“PO  •'>•")  ,  (30) 

f(x)  =  2  TT  2  ® 

p/^ 


f(x)  =  (27r)  \Y(x')e 


-3  f  -  h’)_  , 

^  \  /^/  /vl\  o  fix' 


o  iPn^l  r 

{2tt)  J  a^x^'^e  dx^' 


From  (30)  and  (31)  we  obtain 


■‘Po*l  1/2  "(Pq  '  “  •  P  (-Pq  ’  "  ’  “) 


Ip(*i)  p  “S' 


In  this  case  the  description  obtained  is  only  an  integral  of  the  area  function 
a  thus  yields  very  little  information.  For  =  0  the  left  side 

of  (32)  becomes  the  volume  of  the  target,  but  the  right  side  is  indeterminate. 
However,  the  result  is  of  little  interest  because  of  the  failure  of  the  physical 
optics  approximation  at  low  frequencies. 
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SECTION  V 


15 


where 


a 


P(0,  P2,  P3)  i(P2X2  +  P3X3) 


2  2 
P  +  P 
^2  ^3 


‘“Pz'^Ps 


(38) 


and 


x„  ,  x„)  -  2'^tt  Im  Jx  (0  ,  p^  ,  P3) 


p  X  ') 

^3  3/ 


2  ’  3 


^  C^2^2  ^  ^3”3 -  , 

e  dp^dpg  . 

(39) 


By  inserting  (34)  in  (14)  we  obtain 


f(x)  =  (27r)“^  Jr(x')  (x^'  -  x^)6  (x^'  -  x^ )  ^(x^'  -  x^)  dx' 


=  (27r)-^  (x^-  >  ^2  ’  ^3)  -  ^1’)  • 


(40) 


It  follows  that  (37)  is  satisfied  with 


a 


i^2  ■  ^3)  '  [  >'2  (=*2  ■  ■  ^‘s)  ] 


and 


^  (^2  '  ^3)  '  -  I  ^'**1'  '  -  <'*’^>'^[^2^  (^2.’  ^3)  ■  5'A^2'  ='3')]- 

(42) 

Here  x  =  y  ( x  ,  x  )  and  x  =  y  ( x  ,  x  ')  are  the  functions  that 
1  1'23  1  2^2  3' 

define  the  two  surfaces  enclosing  B  as  illustrated  in  Figure  1. 


16 


Figure  1 .  Functional  Description  of  Body  Geometry 


We  have  tacitly  assumed  that  each  line  parallel  to  the  x^-axis  cuts  the 

2  2 

surface  of  the  target  at  most  twice.  Since  -  y^  =  ^y^  -  y^)  (^yg 
it  follows  that  y  +  y  =  -2/3  /cl  and  y  -  y  =  27r  a  .  These  equations 

^  X  ^  X 

are  easily  solved  to  yield 


=  - 


TT  0!  +  /3 


TT  CL 


a 


CL 


(43) 


Since  a.  and  [3  are  given  by  (41)  ,  (42)  and  (36)  in  terms  of  p(0  ,  p  ,  p  ) 

^  O 

9  /  X 

and  -r —  P  ( 0  ,  p„  ,  p„  )  we  have  obtained  a  complete  determination  of  the 
dp^  2  3 ' 

target  geometry.  We  note  that  the  width  function  is  given  by 
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27r  a 


This  result  agrees  exactly  with  the  result  (22),  (23)  obtained  by  a  different 
method. 

dp 

Our  complete  solution  (43)  requires  a  knowledge  of  p  and  ^  in  the 

"equatorial  plane"  p^^  =  0  ■  Thus  all  frequencies  and  all  aspects  in  this 

0 

plane  are  required.  To  evaluate  (^0  ,  p^  ,  Pg)  we  also  need  aspects 

infinitesimally  near  the  equatorial  plane.  This  result  is  very  closely  related 
to  that  obtained  in  Reference  3  by  an  entirely  different  method.  Since  all 
frequencies  are  required  we  refer  to  the  plane  p^  =  0  as  a  "full-frequency 
plane!'. 
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SECTION  VI 

THE  BAND- LIMITED  PLANE 

2co 

Let  us  recall  that  the  magnitude  p  of  the  vector  p  is  p  =  —  , 
where  cu  is  the  frequency  and  c  is  the  speed  of  light.  Thus  points  near 
the  origin  in  p- space  correspond  to  low  frequencies,  and  distant  points 
correspond  to  high  frequencies.  The  frequency  is  constant  on  spheres 
centered  at  the  origin.  In  most  practical  radar  applications  the  usable  fre¬ 
quency  band  is  limited  by  minimum  and  maximum  values  and  o)^  .  If 
we  set 

2a;  2a; 

m  =  ,  M  =  —  ,  (45) 

c  c 

then  the  usable  portion  of  p- space  corresponds  to  the  annular  region 
m  ^  p  <  M  between  the  two  spheres  illustrated  in  Figure  2.  It  is  im¬ 
portant  to  note  that  the  frequency  limitation  is  due  not  only  to  technological 
factors  but  to  the  failure  of  the  physical  optics  approximation,  on  which  our 
theory  is  based,  at  the  low  frequencies.  Therefore,  in  applying  the  general 
theory  of  Section  III  the  region  D  should  certainly  not  intersect  a  "forbidden 
sphere"  about  the  origin  in  which  the  physical  optics  approximation  fails. 

With  the  exception  of  example  4  of  Section  IV,  all  of  the  cases  examined  so 
far  have  violated  this  criterion,  and  the  exceptional  case  yielded  very  little 
information  about  the  target.  The  "full-frequency  plane"  utilized  in  Section  V 
passes  through  the  origin  and  thus  intersects  the  forbidden  sphere. 

In  this  section  we  will  take  D  to  be  the  disc  consisting  of  that  part 
of  the  plane  p^  =  m  which  lies  within  the  outer  sphere  p  =  M  of  Figure 
2.  The  disc  is  illustrated  by  the  heavy  line  in  Figure  2.  Its  radius  is 
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Figure  2.  Usable  Portions  of  p-Spoce 


Thus  the  region  D  not  only  avoids  the  forbidden  low-frequency  sphere  but 
lies  entirely  within  the  usable  frequency  band.  The  function  K(p)  ,  which 
must  vanish  outside  D  ,  is  taken  to  be 

/  6^p^  -  m)  ,  p  <  M  i 

K(p)  =  }  >  .  (47) 

”  /  2  2  \ 

(  0  ,  p  >  M  ; 
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Then,  from  (12), 


K(x)  =  (27r)  ^  J 


-1 

-dp  ^  "  ^m(^2’  ^3)’ 


2  2 
p 


where 


‘(Va-^Vs),  , 


2  2  2 

P2  -^Ps  'P 


We  note  that  for  M  -*  «> ,  p  approaches  infinity  and  6  approaches 

M 

the  two-dimensional  delta  fimction, 

'’(=‘2 '  ^3)  ”  ^(^2)  H^s)  • 

Thus  6^  is  an  approximation  to  (50).  In  fact,  if  we  set 
(^2  ’  ^3)  "  b(cosi3  ,  sini3)  >  (  ^2  ’  ^3)  "  r(cosQ'  ,  sino;)  ,  /I '  =  /?  - 


we  see  that 


27r  M 


*u(*2’  ’‘3)  '  .W'j’bdbe 


ibr  cosjS 


0  0 


=  (27r )  ^  ^  b  (rb)  db  =  (  27r  r^ )  ^  ^  z  Jq(z)  dz 


^  '^1  " 


/  2  2 
+  X 

2  3 
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Here  and  are  the  Bessel  functions  of  order  zero  and  one,  and  we 
have  used  the  identity 

2 

From  (52)  it  is  easy  to  show  that  6  (0  ,  0)  =  ^  ,  and  the  width  of  the 
large  central  lobe  is  of  order  —  .  Thus  6  will  be  a  good  approximation 
to  the  delta-fionction  (50)  provided 

i  «  L  ,  (54) 

M 

where  L  is  a  typical  target  dimension. 

Assuming  that  (54)  is  satisfied,  we  may  approximate  (48)  by 


imx 

K(x)  =  (2ir)'  e  6(k^)  6(xg). 

Proceeding  with  the  method  of  Section  HI  we  find  from  (13)  that 

_?  _s/o  c  /.  N  e'-  ■  -dp 

f(x)  =  2  TT  \  6  ^p  -  m) - - - 


(irm)  ^  e 


F(m  , 


'^2  ’  ""S, 


where 


(56) 
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F{in,  X2.  Xjj)  j' 


-2-3/2  f  P(“-  P2-P3)*'>  'f’2'  P3) 


J  2  2  2 

2  2  „2  ”  *  P2  *  P3 

P2  ■"P3 


‘(P2*2  P3^3),  , 

X  e 


By  inserting  (55)  in  (14)  we  obtain 


1  r  ■  ^1'^ 

f(x)  =  (27r)  \r(x')e  6  (x^  -  x^')  6{x^  -  x^’)  dx' 


imx 

(27r )  e  ^ 


*1'  ■  '‘2  ■  ’‘3>  P  *'l’ 


We  again  introduce  the  functions  (x^  ,  and  y2C^2  ’  ^3^  illustrated 
in  Figure  1.  It  is  convenient  to  replace  thena  by  functions  a^x^  ,  x^)  and 
b  (X  ,  X  )  defined  by 


yj  =  a-b.y^  =  a  +  b,  a  =  ^(y^+y^).  b  =  1  (y^  -  y^)  =  iw  . 


Here  w  is  again  the  "width  function".  Now,  (58)  yields 


a  +  b 


F(m  ,  x^  ,  Xg)  =  7rme 


1  r  1 

f(x_)  -  2  J  " 


a  -  b 


-ima(x2  ,  x^) 


sin  [mb(x2  ,  x^)] 
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We  see  from  (57)  that  F^m  ,  can  be  obtained  from  measure¬ 

ment  of  p  in  the  disc  D  and  its  image  D'  indicated  by  the  dashed  line  in 
Figure  2*.  Then  the  functions  a  and  b,  which  determine  the  target  completely, 
can  be  obtained  from  (60)  in  a  variety  of  ways.  One  way  is  as  follows.  We 
note  that 

]  F  ]  =  ±  sin  (mb)  (61) 


hence 


^  1^1 


=  ±  b  cos  (mb) 


We  see  that  (62)  consists  of  a  rapidly  fluctuating  "carrier"  modulated  by  the 
slowly  varying  function  b.  Thus  the  (positive)  function  b  is  given  by 

b(x2  .  Xg)  =  envelope  of  1  F  (m  ,  x^  ,  x^  )  1  .  (63) 


Furthermore,  from  (60) 


Hence 


[b  cos  mb  -  ia  sin  mb]  ] 


2  2  2  2 
b  cos  (mb)  +  a  sin  (mb) 


2  •  2  ,  ,  ,  I dF  j 
a  sin  (mb)  =  h; — 
dm 


2  ,  ,  i  rdF 

(^2  '  *3)  '  W 


■  (#)‘| 


In  the  next  section  we  will  attempt  to  eliminate  the  need  to  measure  p  in  D' 
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This  determines  a  (x  ,  \  except  for  an  ambiguity  of  sign  which  can  be 

resolved  by  noting  that 


a 


1  F 

—  arg  - 

m  sin  (mb) 


(67) 


It  is  interesting  to  note  that  the  condition  (54)  for  the  validity  of  our 
approximation  of  6  by  5  can  be  stated  in  the  form 


1  «  Lp 


L 


2 

m 


in  T  ,  2  ^ 

— -  =  LM\I  1  -  cos  6 


2(j0 

=  ML  sin0  =  — -  L  sine  . 
c 


(68) 


Here  co  is  the  upper  frequency  limit  and  6  is  the  semi-angle  of  the  cone 

di 

of  required  aspect  directions  illustrated  in  Figure  2.  We  note  that  0  can 
be  made  arbitrarily  small  provided  co  is  sufficiently  large.  There  is,  of 
course,  no  need  to  take  m  =  2co^/c  to  correspond  to  the  low-frequency 
limit.  We  need  only  require  that  co^  be  greater  than  or  equal  to  this  lower 
limit.  Thus  the  disc  D  will  still  lie  in  the  usable  frequency  band. 
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SECTION  VII 


HALF- TARGET  RECONSTRUCTION 

If  we  think  of  the  p^-axis  of  Figure  2  as  pointing  in  the  direction  of  a 
transmitter  as  seen  from  the  target,  there  are  many  practical  situations  in 
which  one  would  like  to  eliminate  the  need  of  measuring  p  in  the  region  D' 
corresponding  to  aspects  in  the  "back"  of  the  target.  One  might  hope  to 
obtain  a  representation  of  the  "front"  half  of  the  target  without  these  measure¬ 
ments.  The  purpose  of  this  section  is  to  show  that  this  "half- target  recon¬ 
struction"  can,  in  principle,  be  accomplished  by  a  modification  of  our  method. 

For  a  given  target  B  and  a  given  unit  vector  J  ^  we  introduce  the 
symmetric  target  obtained  by  reflecting  the  illmninated  part  (at 
aspect  J^)  of  the  surface  of  B  through  the  origin.  A  typical  case  is  illus¬ 
trated  in  Figure  3.  The  resulting  B^  depends  of  course  on  the  choice  of  the 
origin  o  in  B  . 

It  is  clear  from  (3)  that  if  corresponds  to  B^  and  p  to  B ,  then 

p(p)  «  p^  (p)  for  J  «  .  (69) 

Here  J  =  —  p  =  -  I^  [  see  Equation  (5)]  .  Furthermore,  because  of  the 
symmetry  of  B^  , 

Po(-P)  =  •  (70) 

Using  subscripts  "o"  to  refer  to  B^  ,  we  see  from  (13)  and  (70)  that 

£„(x)  =  [p^(p)  +  (p)]  ■  ?dp  .  (71) 

P 
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Figure  3.  Representation  of  the  Symmetric  Target 


We  now  assume  that  K(p)  vanishes  except  for  ^  -  P  near  . 

Then  from  (60) 

f^(x)  «  ^  [p(p)  +  P*  (P)]  e  ■  ?dp  .  ^^2) 

°  ~  P 

But  from  (14) 

f  (x)  =  V  y  (x')  K(x  -  x')  dx'  =  y  *K  ,  (73) 

and  we  see  that,  by  eliminating  f^(x)  in  (72)  and  (73),  we  obtain  an  equation 
not  for  y  but  for  the  characteristic  function  y  of  B  .  If  we  can  use  this 
equation  to  determine  B  ,  then  we  have  determined  the  illuminated  part 
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of  B  .  The  input  data  in  (72)  requires  the  measurement  of 

P(P)  +  P  (p)  =  2  (Re  p  (p)  (74) 

in  the  region  D  where  K(p)  does  not  vanish.  By  assumption  this  region 
involves  aspects  £  only  near  . 

If,  for  example,  we  apply  the  method  of  this  section  to  that  of  Section 
VI  we  take  _J  =  (1,  0,  0)  to  be  the  unit  vector  in  the  direction  of  the 
p^-axis.  Then  we  merely  replace  p  (-p)  by  p  (p)  iti  (56)  and  (57).  Thus 
we  obtain 


(Re^p  m  ,  p^  ,  p^J 

2  2  2 
m  tp^  H-Pg 


X  e 


+ 


%%  ■ 


(76) 


and  by  using  instead  of  F  in  the  succeeding  equations  of  Section  VI 
we  obtain  not  the  target  B,  but  its  symmetric  replacement  •  We  note 
that  the  input  data  for  (75)  require  the  measurement  of  p  on  the  disc  D 
but  not  on  its  image  D'  . 
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SECTION  vm 


THE  GENERAL  THEORY  FOR  BODIES  OF  REVOLUTION 

In  this  section  we  will  show  that  the  general  theory  developed  in  Section 
III  can  be  considerably  simplified  if  the  target  B  is  axially  symmetric.  In 
effect  we  will  reduce  the  complexity  of  our  formulas  from  three  to  two 
dimensions. 

We  begin  by  introducing  cylindrical  coordinates  in  both  x-space  and 
p- space  defined  by 

X  =  ^x^  ,  r  cos  a  ,  r  sin  oi'j  ,  p  =  ^p^  ,  b  cos  13  ,  b  sin  jS  ^  .  (76) 


Thus 


and 


dx  =  dx^  rdrda  ,  dp  =  dp^  bdbd/3  , 


(77) 


p  •  X  =  +  hr  cos  {13  -  a)  .  (78) 

We  take  the  x^-axis  to  be  the  symmetry  axis  of  the  target.  Then  p(p)  is 
axially  symmetric  (independent  of  ^  )  ,  i.  e. , 


P  =  ,  b)  ,  (79) 

and  the  characteristic  function  is  also  axially  symmetric  (independent  of  a), 
i.  e. , 


(80) 
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If  we  use  the  integral  representation  of  the  Bessel  function 


27r 


J  (z) 
0 


2tt  J 


izcos0 
e  d0  , 


(81) 


and  take  K  =  K(p^  ,  b)  ,  we  find  that  Equation  (12)  becomes 

2  TT  oo 


K 


(x^  ,  r)  -  {2tt)  ^  J  ^ 


ij^p^x^  +  br  cos  (^3  -  q')J 


0  0 


r  r  ^^1^1 

=  (27r)“^  J  dp^  J  bdb  k(p^  ,  b)  e  (br)  . 


(82) 


Since 


(^2 


2  2  2 
^2')  ^  (^3  ~  ^3'}  ~  o!  -  r'  cos  a') 


+  (r  sin  a  -  r'  sin  01') 

2  2 

=  r  +  (r')  -  2rr'  cos  {oi' 


a) , 


(83) 


Equation  (14)  becomes 
f^x^  ,  r^  =  ^dx^'r'dr'do'’'y  ,  r' ^ 

X  K 


X 

-  X  '  , 

-  X  ' 

fx  -  X 

1 

1  ’  ^/ 

2 

2  / 

I 

\  3  3  ;  _ 

2tt 


=  jdx^'jr'dr’y  (  x^'  .  r')  j" 
0 


da' 


0 


X  K 


2  2 

x^  -  X  '  ,  /  r  +  (r')  -  2rr'  cos  a' 

1  1 


(84) 
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Then  from  (84)  we  see  that 


bdb  K I 


iPfX^ 


Jo(br') 


0 

(88) 
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On  the  other  hand,  Equation  (13)  becomes 


t(Xi  ,  r) 


2  \  j  dpj  j  bdb 
0 


‘(P^l 


bdb 


K(Pj  .  b) 

P,^  * 


P  (Pj  ,  b)  +  p*  (-p^ 


X 


iPlX^ 


Jq  (br) 


(89) 


The  application  of  these  formulas  is  much  the  same  as  in  Section  III. 

We  choose  a  domain  D  in  (p^  ,  b) -space  and  a  convenient  function 
K  (p^  ,  b)  which  is  non-zero  only  in  D.  Then  we  determine  Q  from 
Equation  (88)  and  f  in  terms  of  measurable  values  of  p  from  Equation  (89). 
Finally,  we  use  (87)  and  the  fact  that  y  (x^  ,  r)  is  a  characteristic  function 
to  obtain  a  partial  or  complete  determination  of  B.  The  method  may  be  illus¬ 
trated  by  the  following  elementary  examples: 


Example  1: 

D:  all  of  ^p^  ,  b) -space;  K  (p^  ,  b)  =  1  .  (90) 

* 

Since 

OO 

J  bdb  J^(br)  J^(br')  =  ^  d  (r  -  r’)  (91) 

_  0 

*  This  identity  is  equivalent  to  the  inversion  theorem  for  Hankel  transforms. 
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it  follows  that 
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Both  of  the  examples  discussed  here  suffer  from  the  defect  that  p  must 
be  measured  at  forbidden  low  frequencies. 
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